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We calculate bottom-hadron mass splittings with respect to Bd and At using full QCD with 2+1 
, flavors of dynamical Kogut-Susskind sea quarks and domain-wall valence quarks along with a static 

' heavy quark. Our lattices have spatial volume of (2.5 fm)^ with lattice spacing about 0.124 fm and a 

' range of pion masses as low as 291 MeV. Our results are in agreement with experimental observations 

^Nj , and other lattice calculations within our statistical and systematic errors. In particular, we find the 

mass of the Qb to be consistent with the recent CDF measurement. We also predict the mass for 
the as yet unobserved Hj, to be 5955(27) MeV. 
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PACS numbers: 14.20.Mr,12.38.Gc,14.65.Fy 



I. INTRODUCTION 



^ i' Recently, there have been exciting developments in heavy-hadron physics, both theoretically and experimentally. 
O ] Along with discoveries of charmonium states[l|, 0, considerable progress has been made in studying bottom hadrons. 
1-^' While the B factories, such as Belle and Babar, have investigated the bottom mcsons[l|, recent experiments at 
Fermilab have reported the discoveries of a few bottom baryon states. In the summer of 2007, CDF[| reported the 
fSJ i first observation of the heavy baryons St and SJ, and then both D0[1] and CDF0] observed the bottom baryon 
J> ' S^, breaking a long period of silence following the observation of the A^ in 1991. The bottom baryon spectrum 
. has become somewhat controversial due to recent results from DO and CDF. Last summer, DOQ reported a first 
' observation of the doubly strange bottom baryon fl^ at 6.165(10) (13) GeV. However, a recent CDF work puts the 
fi^ mass at 6.0544(68)(9) GeV, a difference of 111(12)(14) MeV. With a discrepancy of 6.2 standard deviations, it 
appears that the two collaborations cannot both be observing the Qb- A theoretical understanding of the bottom 
baryon spectroscopy from first-principles QCD is crucial and can help disentangle such discrepancies in experiment. 
Furthermore, it is anticipated that in the upcoming dedicated bottom physics experiment LHCb at CERN. there will 
I be many more discoveries in the bottom-hadron spectrum. Combining such experiments with improved theoretical 
understanding, our knowledge about these states will in the near future be significantly enhanced. 

Lattice QCD has been successfully computing hadron masses for the last few decades with increasing precision and 
scope. However, the study of heavy quarks requires special techniques, such as the introduction of an action that 
5^ ■ will minimize 0{ma) errors due to the discretization of space. Such errors are more prominent on coarse lattices due 
1 to the heavy quark masses involved. Many existing bottom spectroscopy lattice calculations have been performed 
in the meson sector. Ref. calculated the bottomonium and bottom-light meson spectra using the relativistic 
heavy-quark action Q with fermion action parameters determined through nonperturbative tuning to eliminate 
systematic uncertainties; their ensembles used 2-1-1 flavors of domain- wall fcrmions (DWF)(Tol. [Tl]| with lightest pion 
mass 275 MeV. Ref. [l2| also used DWF gauge ensembles but with NRQCD[i^ to calculate bottom-quark quantities. 



A more impressive work|l^ done by the TriLat Collaboration used anisotropic 2-fiavor Wilson-type lattices with 
renormalized anisotropy as high as 6 and operators projected into irreducible representations of the lattice cubic 
group. These techniques allowed them to obtain very clean signals including multiple static-light excited states. A 
more complete review of previous studies can be found in the lattice review talks of Refs. [l^, QQ, [13, and 
references therein. 

Unlike bottomonium and B mesons, bottom baryons have not received as much attention from lattice QCD. Some 
pioneering worksfU [2^ were done using extrapolations with light-fermion actions or the NRQCD action in the 
quenched approximation, where the fermion loop degrees of freedom are absent. Since it is difficult to estimate the 
systematic error due to the quenched approximation, high-precision calculations cannot be achieved. Recently, more 
dynamical ensembles have become available due to the increase of the computer resources available for numerical 
research, and since the recent discovery of the double-5 baryon, more lattice calculations have emerged. Ref. 
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calculated single- and doublc-6 baryons with NRQCD action for the bottom quark on isotropic 2-|-l-flavor clover 
ensembles, having lightest pion mass around 600 MeV. Ref. [2^ went to a lighter pion mass (275 MeV) using 2+1 
flavors of domain-wall fermions but used static-quark action to simulate the bottom quarks. Since the static-quark 
action tends to result in noisy signals, one needs high statistics to yield numbers comparable with experimental results. 
A selection of single-6 baryons were calculated using static-quark action on 2-flavor chirally improved lattice Dirac 
operator at pion masses as light as 350 MeV in Ref. An ongoing work using staggered and Fermilab fermion 

actions [131 on MILC lattices with multiple lattice spacings was presented in Ref. [28| . 

In this work, we report the mass-splittings of the bottom-hadron spectrum using a static-quark action to simulate 
the bottom quark and using domain- wall fermions for the light valence quarks. We use an extensive set of gauge 
ensembles of 2-|-l-flavor staggered-fermion lattices with a range of quark masses resulting in pion masses as light 
as 290 MeV; the number of available configurations for these ensembles allows us to achieve high statistics. These 
ensembles have lattice spacing a = 0.124 fm, and the lattice volume is about (2.5 fm)"^. It is anticipated that the 
physical bottom spectrum resembles the spectrum of hadrons with a single static (infinitely massive) quark, where 
corrections are suppressed by powers of AqcD/^if,, which can be included systematically in heavy-quark effective field 
theory. 

II. LATTICE SETUP 

The gauge configurations used in our work were generated by the MILC collaboration using the one-loop tadpole- 
improved gauge action[29|, where both 0{a?) and Oig^o?) errors are removed. For the fermions in the sea, the asqtad 
improved Kogut-Susskind action^, jsil, 32, 33, 34, 33| is used. This action is the Naik actionfsBl (O(a^)-improved 
Kogut-Susskind action) with the one- link terms smeared such that couplings to gluons with any momentum component 
equal to IT /a are set to zero. 

For the valence-sector light quarks (up, down and strange), we use the five-dimensional Shamir domain- wall fermion 
action pH Wf\. This action introduces a fifth dimension of extent L5 and a mass parameter M5. In our case, we used 
L5 = 16 and M5 = —1.7. The physical quark fields q(x,t) reside on the four-dimensional boundaries of the fifth 
dimension. The left and right chiral components are exponentially separated onto the opposing boundaries, resulting 
an action with approximate chiral symmetry at finite lattice spacing. The bare-quark mass parameter 
introduced as a direct couphng of the boundary chiral components. Hypercubic-smeared (HYP-smeared) j38l[39[|40l . l4l| 
gauge links were used in the domain- wall fermion action to improve chiral symmetry. 

Because the valence-quark and sea-quark actions are different, the calculation is inherently partially quenched. 
That is, the calculation violates unitarity. Unlike conventional partially quenched calculations, which become unitary 
when the valence-quark mass is tuned to the sea-quark mass, unitarity cannot be restored by tuning at nonzero 
lattice spacing. The next-best option is to tune the valence-quark mass in such a way that the resulting pions have 
the same mass as those made of the sea Kogut-Susskind fermions. In this case unitarity should be restored in the 
continuum hmit, where the nf = 2 staggered action has an SU{8)l ® SU{8)ii U{l)v chiral symmetry due to 
the fourfold taste degeneracy of each flavor, and each pion has 15 additional degenerate partners. At flnite lattice 
spacing this symmetry is broken, and the taste multiplets are no longer degenerate but have splittings that are 
0(Q;^a^)[33, HH, m, [3^, E^l. The domain- wall fermion mass is tuned to give valence pions that match the Goldstone 
Kogut-Susskind pion. This is the only Goldstone boson that becomes massless in the chiral limit at nonzero lattice 
spacing. This choice gives pions that are as light as possible, resulting in better convergence of the chiral perturbation 
theory (XPT) needed to extrapolate the lattice results to the physical quark mass. This tuning was performed and 
used by LHPC collaborationlillil,!!!,!!!,!!^,!!! as weU. 

The gauge ensembles used were generated for five pion masses down to 290 MeV with 20"^ x 64 lattice volume, 
corresponding to a physical volume V = (2.5 fm)^. The lattice spacing is a = 0.12406 fm, as quoted in Ref. [i^, 
which was determined from T spectroscopy. Further information about our lattice actions, parameters and tests of 
our computational setup can be found in Ref. [soj . The details concerning the ensembles used in this calculation can 
be found in Table [D 

A. Static-Quark Action 

Starting from the Wilson fermion action on the lattice, we can derive the quark propagator at the static limit, 
where the quark mass is infinite. In this limit the quark mass is removed, leaving only the forward hopping term in 
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Label 


IVU (GeV) A^conf 




m007 


0.2938(1.1) 


478 


7 


m007' 


0.2938(11) 


429 


24 


mOlO 


0.3566(08) 


657 


6 


m020 


0.4969(06) 


484 


7 


m030 


0.5982(08) 


564 


5 



TABLE I: The gauge-ensemble parameters, including the pion mass Af,r in GeV, the number of gauge configurations A'conf and 
numbers of sources A'src used in our calculation. The m007' configurations differ from m007 only in the number of sources used. 

the action: 

GQ(x,i;io) = ^^ n t^4(x,t')- (1) 

The links U4 that enter this propagator can be any set of gauge covariant paths that connect neighboring points in 
time. Different choices of paths result in actions that arc equivalent up to discretization errors. All of them have 
the same continuum limit. In our calculation, we tested two choices for the temporal links entering the static quark 
propagator. We used the simple temporal gauge links of the underlying gauge configuration. In addition, we used 
HYP-smcarcd links with the same smearing parameters as those used in the domain- wall action (38l|. We observed 
that the HYP-smeared links gave significantly better signals for our correlators; hence, we only report results from 
that choice. 



B. Static-Light Correlation Functions 

At the static-quark limit, the heavy-quark propagator is reduced to a Wilson line. We can construct the static- 
light hadron correlation functions by contracting the static-quark line with light-quark propagators to create gauge 
invariant correlation functions. The static-light baryon correlation functions have the form 

Gr(x,t; xo,io) = t)rgj,(x, t)eca6^"''''(x, xq, io)g/'(^o, to)rg/ (xq, to)ec'a'&'), (2) 

and the static-light meson correlator is 

M{x,t;xo,to) = {q'}{x,t)P'"''{x,t;xo,to)q'}'{xo,to)), (3) 

where / G {u, d, s} is the light-quark flavor index and P'^'^ is the Wilson line connecting the source and the sink which 
are separated by time t. The spin matrix F is either C75 or 6*7^ for scalar or vector diquarks, respectively. 

The quantum numbers (spin J, isospin / and strangeness S) of the static-light hadronic correlators we construct 
are listed in Table [Til Since the static quark carries no spin or light-flavor quantum numbers, the numbers in the table 
correspond to the quantum numbers of the diquark. Since the gamma matrix insertions in the baryon correlators 
are symmetric (C7^) or anti-symmetric (C75), the corresponding flavor wavefunctions arc symmetric (sextet) or 
antisymmetric (antitriplct), respectively. This flavor symmetry structure is automatically enforced without explicit 
symmetrization or antisymmetrization of the flavor wavefunctions. In the static limit the hyperfine splittings are 
exactly zero, so we cannot compute such spin splittings. For the physical bottom quark the hyperfine splittings are 
less than one percent. When comparing to experiment, we identify our states with the lowest-spin physical particles. 

III. NUMERICAL RESULTS 
A. Effective Mass Plots 



The hadron spectrum can be calculated in lattice QCD using two-point hadronic correlation functions 

C{to,t) = (O|0(t)0t(io)|o), 



(4) 
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State 


J I s 


Bi 


1/2 1/2 


Bs 


1/2 -1 


Ab 





Si, 


1 1 




1/2 -1 




1 1/2 -1 


ilb 


1 -2 



TABLE II: The (diquark) quantum numbers of the states studied in this work. The names used are motivated by the flavor 
content of the states. 



using the creation and annihilation operators and O described in Sec. [TTl In order to extract the masses of the 
states of interest, we take momentum (= to extract the masses) and spin projections at the source and sink. By 
inserting a complete set of hadronic eigenstates of the Hamiltonian between the creation and annihilation operators 
we get 



(5) 



where An = |(0|O|n)p is the overlap factor between the n'^ eigenstate and the state created by the operator, and Af„ 
is the mass (or energy if considering a multi-particle state) for the n'^ eigenstate. In addition to differences in the 
flavor and spin structure of the operators, we also employ smearing on the fermion fields to improve spatial overlap 
with the hadronic wavefunction. We present results using smeared operators at both source and sink (SS), and also 
using smearing only at the source with a point sink (SP). 

At large time t, the states with higher masses become exponentially suppressed, and Eq. [5] becomes dominated by 
the ground state with mass Mq. Thus, we expect the correlator at large time to proportional to e~^^°*. A simple 
manipulation allows us to cancel the overlap factors and extract the "effective mass" at each t: 



Meff (t + 1/2) = log{Cit)/C(t + 1)). 



(6) 



For sufficiently large lattice time extent and sufficiently small noise, Mgg{t) will asymptotically approach Mq at large 
time. In practice, there will exist a window in which the effective mass is approximately flat, called the plateau, 
between the early times which are contaminated by excited states and later times which are dominated by statistical 
noise. By applying more sophisticated manipulations to the correlators, it is possible to address excited states in the 
effective mass analytically using single or multiple correlators as input; for more details, see Refs. [sil. Is^. Issl. Is^. Issl. IsgI . 
Since we are only interested in the ground state in this work, Eq. [S]is sufficient to crosscheck the fitted masses. 



0.7 



0.6 



5 



0.5 



0.4 



I ' I ' I ' I ' I ' I ' I ' I ' I ' I ' I ' I ' I ' I ' I ' I 



^ @ © jgC at gj 



1 5 5 J 



I . I . I . I . I . I . I . I . I . I . I . I . I . I . I . I 



10 12 14 16 



FIG. 1: Effective mass of the lightest Bg. Note the oscillation at small time. Both axes are in lattice units. 
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B. Fitting Metliodology 

111 the (Euclidean) continuum, the two-point correlators follow the form Eq. [5l We can truncate this expansion at 
the first term as long as we apply such a single-exponential fit only to the plateau region of each two-point correlator. 
However, since we are using the DWF action in the light-valence sector, an additional complication arises from the 
non-locality of the DWF action in four dimensions. This non-locality manifests as oscillations in the effective mass 
close to the source, as seen in Fig. [T] in the range 3 < t < 7. A phenomenological form for fitting such data was 
proposed in Ref. [46|, employing an oscillating contribution describing the non-local artifacts: 

Cit) = • • • + Ao,e cos(7rOe-^^-(*-*°), (7) 

where the dots indicate that this is a correction to the continuum form. Since we find empirically that Mqsc is quite 
large compared to Mq, at large t the oscillations are suppressed and the ground-state mass becomes dominant. 
The mass Mh of a hadron with one heavy quark can be expanded in powers of the heavy quark mass: 

Mh = ZqiTiQ + A + — + ••■ , (8) 

mQ 

where mg is the mass of the heavy quark and Zq, A, and c are arbitrary coefficients. In the static limit the leading 
term proportional to the heavy-quark mass is eliminated while the sub- leading corrections proportional to l/mg, 
although calculable in the heavy-quark effective theory (IIQET)[53|, are ignored. For this reason, in the static limit 
only mass splittings between hadrons with the same number of heavy quarks can be easily computed. For two hadrons 
Bi and B2 the leading term in the ratio of their two-point correlators is 

{Cb, {t)/CB. it)) = A[,e-^^^"(*-*") + a;,, cos(^t)e-^<=o(*-*°) + • • • , (9) 

where the dots indicate the effects of the excited-state terms which are suppressed exponentially at large time. The 
lightest meson and lightest baryon (A;,) are natural candidates for the denominator of correlator ratios. In two 
special cases, we try other combinations in order to make more direct comparisons to previous lattice results. 

We will compare two kinds of fits: a simple single-exponential and a single-exponential plus the oscillating term. 
The use of the oscillating term may allow us to extend the fitting range to earlier times, but it comes at the cost of 
introducing twice as many parameters. We determine the appropriate fit range for our fits using a two-step heuristic: 
First, calculate the Q value associated with constant fits to the effective mass over all possible ranges. Then, apply 
the appropriate fit (single-exponential or oscillating) to that range. If extending the range to earlier or later times 
does not diminish the Q value by more than a factor of 2, extend the range. Continue to extend to the maximum 
length possible. This heuristic provides a compromise between optimizing the quality of the fit and ensuring that the 
fit ranges include as much of the data as possible. 

Once the appropriate fit range for each correlator ratio is determined, we apply a Levenberg-Marquardt algorithm 
to minimize for each fit form. The procedure is carried out under a single-elimination jackknife to estimate the 
error in the parameter determinations. We compare the results of fitting the SP two-point correlators using either 
single-exponential or oscillating fit form in Table Hill Since we see good agreement between these two sets of correlators 
and no uniform improvement by using the more complicated oscillating fit, we adopt the simpler one-state exponential 
and apply it simultaneously to both SS and SP correlators; the results are shown in Table HVl and a selected fits and 
its corresponding effective mass plots arc also shown in Fig. [2l 

C. Chiral Extrapolations 

Having extracted the static-light hadron splitting spectrum for a number of unphysically heavy light-quark masses, 
wc can extrapolate to determine the spectrum at the true light-quark mass. (The strange-quark mass already falls in 
the range accessible to lattice calculations; we use a single value for it and will assign a systematic error associated 
with the mismatch between our adopted mass and the physical strange-quark mass.) For sufficiently light masses, 
it should be possible to employ chiral perturbation theory (XPT) as an extrapolation form. The subject of XPT 
in the presence of heavy quarks (HXPT) has been addressed in numerous works [58,, ^62, ilS,, iG^l, where 

Goldstone bosons and photons are the elemental degrees of freedom. HXPT incorporates both the approximate 
chiral symmetry as well as the heavy-quark symmetry of QCD and have non-analytic terms in addition to the usual 
polynomial terms in pion and kaon masses. Such forms could in principle be used to extrapolate lattice data and 
determine the parameters induced in the theory. In addition, such chiral effective theories can be modified to include 
the mixed-action corrections in the vacuum [65l. [66l. [67j . 
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0.007 0.01 0.02 0.03 






0.007 0.01 0.02 0.03 


Bs — Bd 


0.039(01) 0.043(03) 0.026(01) 0.017(01) 




Bs 


-Bd 


0.039(02) 0.043(03) 0.026(01) 0.017(01) 




0.253(08) 0.298(08) 0.300(12) 0.332(10) 




Ab 


-Bd 


0.253(08) 0.299(08) 0.282(18) 0.332(10) 


Sfc — Bd 


0.329(06) 0.355(07) 0.332(11) 0.365(08) 






-Bd 


0.329(06) 0.355(07) 0.318(15) 0.365(09) 


Sb — Bd 


0.438(09) 0.403(21) 0.411(13) 0.437(08) 






-Bd 


0.438(09) 0.403(21) 0.391(20) 0.437(08) 


Sf, — Bd 


0.456(08) 0.466(08) 0.425(15) 0.455(07) 






-Bd 


0.456(08) 0.466(08) 0.426(15) 0.455(07) 


— Bd 


0.484(07) 0.505(06) 0.484(05) 0.473(06) 






-Bd 


0.485(09) 0.505(06) 0.477(07) 0.473(06) 


Ab-Bs 


0.214(08) 0.255(10) 0.256(16) 0.313(10) 




Ab 


-Bs 


0.214(08) 0.243(13) 0.256(16) 0.313(10) 


Eb - Ab 


0.099(04) 0.068(05) 0.045(01) 0.033(02) 




Ht 


-Ab 


0.099(04) 0.070(04) 0.045(02) 0.033(02) 


- Ai, 


0.190(17) 0.166(09) 0.106(11) 0.120(06) 




Si, 


-At 


0.211(07) 0.177(07) 0.106(11) 0.126(04) 


H'i, - A, 


0.218(14) 0.192(08) 0.138(10) 0.144(08) 






-Ab 


0.234(06) 0.201(06) 0.138(10) 0.141(04) 


fifc - Ab 


0.256(13) 0.224(08) 0.172(09) 0.162(08) 






-Ab 


0.266(06) 0.231(06) 0.182(04) 0.157(04) 



TABLE III: The splittings on ensembles with various pion masses from one-exponential (top) and oscillating fits (bottom) to 
smeared-point correlators. Numbers are given in lattice units. 





0.007 0.01 0.02 0.03 


Bs 


-Bd 


0.040(01) 0.046(03) 0.027(01) 0.018(01) 


Ab 


-Bd 


0.232(06) 0.282(08) 0.279(06) 0.332(10) 




-Bd 


0.318(05) 0.342(05) 0.333(05) 0.352(05) 


Si, 


-Bd 


0.420(09) 0.390(16) 0.386(11) 0.430(08) 




-Bd 


0.441(06) 0.452(08) 0.424(08) 0.448(07) 


fif, 


-Bd 


0.477(07) 0.496(06) 0.466(06) 0.467(06) 


Ab 


-Bs 


0.190(06) 0.230(07) 0.252(06) 0.315(10) 




-Ab 


0.101(04) 0.073(04) 0.052(02) 0.034(01) 




-Ab 


0.177(13) 0.155(06) 0.121(07) 0.109(10) 




-Ab 


0.213(09) 0.175(11) 0.142(08) 0.128(09) 


fit 


-Ab 


0.258(08) 0.224(08) 0.179(07) 0.149(08) 



TABLE IV: The splittings on ensembles with various pion masses from one-state exponential fits to smeared-point and smeared- 
smeared correlators simultaneously. Numbers are given in lattice units. 



However, given the small number of data points available and the large number of parameters in a typical XPT form, 
a comprehensive chiral extrapolation is not possible. Therefore, we adopt only the leading-order form linear in Af^. 
Since HXPTjssj suggests that the next-order non-analytic term can be approximated by a term proportional to A'l^, 
we will also test that form to estimate the systematic error due to uncertainty in the light-quark mass extrapolation. 

Another consideration prompted by the gauge ensembles we have selected is the variation of the lattice spacing 
between ensembles. Each ensemble uses different light sea-quark mass and different gauge coupling, resulting in 
variation of the lattice scale. In order to avoid these ambiguities in the scale determination, we have chosen to 
use dimensionless quantities in the chiral extrapolations. Three different scales, which may be sensitive to different 
physics, are taken in ratios with the mass splittings and pion mass: the pion decay constant /^r, the baryon mass Mn 
and the Sommer scale of the static-quark potential ri; their values for these ensembles can be found in Refs. [iol-lsoj. 
See Fig. [3] for examples of the extrapolations. 

We summarize the extrapolated mass splittings in Table |V] and Figs. 21 We find that there is good consistency 
between the reference scales f-^ and Mn, but that ri is slightly away (< 2a) from the other two. This discrepancy 
should probably be interpreted as a systematic error due to the nonzero lattice spacing. There is also a discretization 
systematic error due to the single lattice spacing used. One needs to calculate at least two different lattice spacing 
to be able to extrapolate to the continuum limit. The differences in using methods for setting the lattice spacing 
are likely to be caused by the lack of such an extrapolation. The discretization systematic can be also estimated by 
examining the quark action. For example, in our light-valence sector, the action has 0(a^ Aq^p) error, which is about 
4%, smaller than the discrepancies we observed here. Therefore, we are not underestimating such a systematic effect 
here. Since Mq has smaller dependence on the light-quark mass and is a well measured experimental state, we will 
use the extrapolated mass splittings normalized by Afo as main estimates of the central values and statistical errors. 
The discrepancies with other reference scales, we use for estimation of discretization systematic errors. 
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FIG. 2: A selection of effective mass splittings showing simultaneous fits to smeared-smeared (dark band) and smeared-point 
(lighter (red) band) correlators. The (blue) triangles and (red) circles in each panel are the smeared-smeared and smeared-point 
effective mass splittings. Both axes are in lattice units. 







ri 


Mn 




-Ba 


71.0(2.3)(1.5) 


78.4(2.4)(1.3) 


71.2(2.2)(1.2) 


At 


-Bi 


326(11)(11) 


372(12)(9.1) 


340(11)(8.1) 


Hi, 


-Ba 


470.8(7.9)(8.1) 


532.6(8.5)(7.0) 484.7(7.7) (6.2) 




~Bd 


600(15)(16) 


676(17)(13) 


615(15)(12) 




-Ba 


656(10)(12) 


739(11)(10) 


672(10)(9.1) 




-Bd 


731(10)(12) 


824(11)(10) 


749.2(9.8) (9.0) 


At 


-Bs 


248(11)(11) 


286(12)(9.5) 


261(10)(8.5) 




-Ab 


157.5(5.2)(4.3) 173.0(5.7)(3.9) 157.2(5.2)(3.4) 




-At 


272(14)(17) 


301(14)(14) 


274(13)(13) 




-Afc 


333(15)(14) 


369(16)(12) 


335(15)(10) 




-At 


411(12)(12) 


456(13)(11) 


414(12)(9.5) 




-Ei, 


61.9(2.7)(2.6) 


69.0(2.9)(2.2) 


62.6(2.6)(2.0) 




1—/ 


81.1(2.8)(4.3) 


89.8(3.0)(3.9) 


81.5(2.7)(3.4) 



TABLE V: The extrapolated splittings using leading-order HXPT with the lattice scale set according to the pion decay constant 
f-n, the Sommer scale ri and the baryon mass Mn- The numbers are given in MeV; first error is statistical and the second 
is a systematic associated with the mass extrapolation. 



We compare our results with previously published lattice calculations using 2-1-1 dynamical flavors. For example, 



Detmold et al.|25j contains all the mass splittings calculated in this work, and they use a similar lattice spacing. 



We find good agreement with their numbers, even though we use a different fermion discretization in the sea sector 
(staggered versus DWF). 

Unfortunately, other 2-|-l-flavor worksjUll^ calculated fewer mass splittings, so we cannot have a comprehensive 
comparison. Lewis et al.jlJj used the NRQCD action to simulate their bottom quarks on CP-PACS clover lattices; 
Na et al.[2^ calculated using the same MILC lattices used in this work but with Fermilab heavy-quark action and 
staggered light-quark action; they also included another two lattice spacings. All four calculations include the splittings 
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FIG. 3: Selected mass-splitting extrapolations with reference scales (top), ri (middle) and Mq (bottom) 



Sfe — Af, and Ef, — At,; we calculate the additional splittings SJ^ — T,b and fib — in order to make direct comparisons. 
The results are shown in Fig. [6l We again see good agreement amongst all lattice calculations, and mild scaling in 
the Na ct al. results. Agreement with experimental values as given by the PDG is fairly good. The discrepancy in 
the Sb — Af, splitting may be due to discretization effects, as suggested by the trend visible in the Na et al. results. 

We summarize our consensus results for the splittings of the bottom hadron spectrum in Table IVTl We adopt the 
leading-order HXPT extrapolation in terms of the AIq scale as our central values and explicitly give three sources 
of error: statistical, extrapolation and scale. The extrapolation systematic error is estimated by the discrepancy 
between the leading-order extrapolation and the result using a form including the next-order term. These errors 
are generally about the same size as the statistical error. The error due to scale setting and discretization is estimated 
using the spread amongst the three scale-setting methods: /,r, fi and Mq. This systematic is quite large, up to four 
times the statistical error, and could be resolved by repeating this calculation on finer lattices. The finite-volume 
correction are negligible for 2.5 fm box, since no such correction is observed in the light-hadron masses (50| either. 
The remaining systematics, such as the effects of the Aqcb /nib corrections are omitted, since they are substantially 
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FIG. 4: Comparison of mass-splitting extrapolated values using different reference scales. Splittings with respect to Bd are on 
the left; those with respect to A;, are on the right. The solid (red) bars indicate the experimental values given in the PDG, 
where available. For the Qb, we show both the DO resultQ] (upper-right, purple) and the CDF result[6^ (lower-left, magenta) 



smaller than our main systematics discussed in the text (less than one percent). 

These results are in fairly good agreement with experimental results. The largest discrepancies appear for the Sj 
although it is only about one sigma away, given the estimated systematics. We give both the D0@ and new CDF[6j 
results for the fib mass. Our result (like most theory predictions) is more consistent with the lower CDF value. 



IV. SUMMARY AND OUTLOOK 



In this work, we calculated the mass splittings of bottom baryons using 2-|-l-flavor staggered-fermion lattice gauge 
ensembles, domain- wall fermion light valence quarks and static-quark action for the bottom quarks. This is a relatively 
high-statistics calculation resulting statistical error bars smaller than other recent lattice calculations of baryon mass 
splittings in the static limit. Our dominant errors are systematic. The most important systematic error is scale 
setting. We have used three options for converting to physical units and the spread of our final results is an indication 
of such errors. Such discrepancy is likely to be caused by a nonzero lattice-spacing systematic error in our calculation. 
A rough estimate of the lattice-spacing systematic error is 4% of the quoted numbers, assuming ©(a^Ag^^) scaling 
violations caused by the light quarks; however, more generous systematic errors are assigned based on the discrepancy 
from setting the scale using different quantities. Extrapolations to the physical light-quark mass also introduce 
errors that should be of the same order of magnitude as typical lattice calculations with the same parameters as 
ours. Finally, Aqcd/"t-& corrections (estimated to be 1%) have to be included for making direct comparison to 
experiment. Nonetheless, all the above systematic errors are small, and comparisons with experiment can provide 
useful conclusions. Our results with our estimates of systematic errors are summarized in Table IVll 

Our results agree well with experiment in the cases of known mass splittings. In the case of f^h our calculation is 
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FIG. 5; A comparison of our mass splittings (using Mn reference scale values) with those of Detmold et al.[231. The error 
estimates for both works are statistical only. The solid (red) bars indicate the experimental values given in the PDG, where 
available. For the Qt, we show both the DO result^] (upper-right, purple) and the CDF resultjg^ (lower-left, magenta) 



in agreement with the recent CDF result and several standard deviations away from the DO result. Such conclusion 
holds for all other lattice calculations of the ili, mass. Our results for the SJ, mass splittings are a prediction since S'j, 
has not yet been observed. Using the splitting from the Af, and summing with the experimental value of that baryon, 
we predict the mass of the SJ, to be 5955(27) MeV, where the error given combines in quadrature all our errors. 

In the future we plan to extend these calculations using the anisotropic gauge ensembles generated by the Hadron 
Spectrum Collaboration (HSC)[6i,[70l. In this case, taking advantage of the anisotropy, a relativistic improved action 
for the bottom quark may be used, eliminating the need for resorting to the static approximation. As a result, 
statistical errors are expected to be dramatically improved, and high precision can be achieved, as demonstrated in 
a recent workfs^l by NPLQCD Collaboration. In this framework, it will also be possible to compute the spin-3/2 
to spin-1/2 mass splittings; such effects are inaccessible in the static limit. In addition we expect to have access to 
lighter pion masses and possibly a second lattice spacing, addressing all sources of systematic error in this calculation. 
A recent paper by HSC[7l| demonstrates a new technique, optimized for complicated operators, which may further 
improve the calculations and in particular allow for the extraction of the excited-state spectrum of the bottom-baryon 
sector. 
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Splitting(Stat.) (Extrap.) (Scale) 


Experiment 






71.2(2.2)(1.2)(4.2) 


87.1(0.6) 


Ab 


-Bd 


340(11)(8.1)(24) 


341.0(1.6) 




-Bd 


484.7(7.7)(6.2)(32) 


513(3) 




-Bd 


615(15)(12)(40) 


554(3) 


—6 


-Bd 


672(10)(9.1)(44) 


- 


fib 


-Bd 


749.2(9.8)(9.0)(49) 


786(7)/886(16) 


Afc 


- B,, 


261(10)(8.5)(19) 


253.9(1.7) 




-At 


157.2(5.2)(3.4)(9.0) 


172(3) 


St 


-At 


274(13)(13)(17) 


213(3) 


—6 


-At 


335(15)(10)(20) 






-At 


414(12)(9.5)(25) 


445(7)/545(16) 




-St 


62.6(2.6)(2.0)(3.9) 




Ob 


— — t 


81.5(2.7)(3.4)(4.9) 





TABLE VI: Our consensus heavy-hadron splittings in units of MeV. The first error is statistical, the second is a systematic 
associated with the mass extrapolation and the third is a systematic associated with the setting of the lattice scale. 
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FIG. 6: Comparison of mass splittings with all available 2+1-flavor lattice calculations of bottom baryons. The square (blue) 
points are the points extrapolated using the f2-mass reference scale; the solid error bars indicate the statistical error, and the 
dashed bars indicate the total errors (including the estimated systematic ones). The solid (red) bars indicate the experimental 
values given in the PDG, where available. For the Ot, we show both the DO resultQ] (upper-right, purple) and the CDF 
result [6a| (lower-left, magenta). 
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